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We examine whether it is possible to describe simultaneously the experimental 
data from RHIC and LHC on nuclear suppression of light hadrons and non-photonic 
single electrons in the pQCD picture of parton energy loss. We perform calculations 
accounting for both radiative and collisional energy loss. We show that once the 
coupling constant is fixed from comparison with data on the nuclear modification 
factor for light hadrons it gives a satisfactory agreement with data on the electron 
Q , Raa and azimuthal anisotropy V2- Our results show that the collisional mechanism 

,-J^ I is only of marginal significance in nuclear suppression of single electrons except for 

the bottom contribution at momenta < 6 -^ 8 GeV. 

> 

^ ■ 1 Introduction 

^ I A remarkable result of the experiments on AA-coUisions at RHIC [1, 2, 3] and LHC [4] 

^ ' is the observation that nuclear suppression of high-p^ single electrons from semi-leptonic 

ff^ ■ decays of heavy mesons is almost as strong as that of pions. It is believed that nuclear 

suppression of high-p^- particles (jet quenching) is due to radiative [5, 6, 7, 8, 9, 10] and 
collisional [11] parton energy loss in the hot quark-gluon plasma (QGP) produced in the 
^ I initial stage of AA-collisions. The observed strong suppression of single electrons indicates 

that heavy quarks are quenched approximately as light ones. This seemed to be somewhat 
puzzling in light of the prediction of the dead cone reduction of the heavy quark radiative 
energy loss [12]. In recent years nuclear suppression of non-photonic single electrons 
due to radiative and collisional heavy quark energy loss in pQCD received considerable 
theoretical attention [13, 14, 15, 16, 17, 18, 19, 20]. However, no clear consensus has 
emerged on the physical picture of the effect and whether the pQCD can explain it. 

Comparison of the radiative energy loss (calculated within the light-cone path integral 
(LCPI) approach [7, 21]) and the collisional energy loss calculated with the same a^ and 
the Debye screening mass shows [22] that for quark energy E >^ 5 GeV the collisional 
energy loss is relatively small for light quarks and c-quark, and for 6-quark the collisional 
energy loss becomes important only at -E < 10 GeV. This shows that the collisional 
energy loss should be of only marginal significance in suppression of single electrons at 
Pt ^ 5 GeV (since the electron spectrum is controlled by heavy quark production at 
approximately twice the electron momentum). This is also supported from computations 



^ 



of the nuclear modification factor Raa of single electrons for purely collisional mechanism 
in [19, 20], where it was found that to fit the data the cross sections of the 2 — )■ 2 processes 
should be enhanced by a factor ~ 4. Thus, it is natural to expect that, if the pQCD is valid 
for parton energy loss for RHIC and LHC conditions, the nuclear suppression of single 
electrons should be described in a picture where the radiative mechanism dominates. 

The purpose of the present work is to analyze the available RHIC and LHC data on 
the electron suppression, and to examine whether the observed nuclear suppression of 
single electrons and light hadrons can be described simultaneously in the pQCD. In our 
study we use the LCPI approach [7] to induced gluon emission. The advantage of this 
formalism is that it treats accurately the mass and finite-size effects, and is valid beyond 
the soft gluon approximation (used in previous analyses [13, 15, 16, 18]). Calculations 
beyond the soft gluon approximation are especially desirable for c-quark. Indeed, in the 
LCPI formalism [7] the induced gluon a:-spectrum {x is the gluon fractional momentum) 
is expressed through the solution of a two-dimensional Schrodinger equation in which the 
longitudinal coordinate z (along the fast parton momentum) plays the role of time. And 
the Hamiltonian (see Appendix) depends on the parton masses only through the term 
1/Lf = [rrigX^ + m^(l — a;)^]/2a;(l —x)E, where E is the initial quark energy, rriq^g are the 
quasiparticle parton masses. The quark mass becomes important when m'^x'^ > m^(l— x)^. 
Taking rrig ~ 400 MeV [23], we see that it occurs at a; > 0.3 for c-quark and a; > 0.1 for 
6-quark, and below these values the sensitivity to the quark mass should vanish quickly. 
Accurate numerical calculations [24] corroborate these qualitative estimates. Note that 
the results of [24] show that at energy ~ 10 ^ 20 GeV for a finite-size plasma the quark 
mass suppression of the radiative energy loss may be considerably weaker than predicted 
in the dead cone model [12], and at higher energies the radiative energy loss may even be 
enhanced for heavy quarks. Therefore the observed strong suppression of single electrons 
does not seem to be very strange. 

We treat the effect of parton energy loss on the single electron yield within the scheme 
developed previously for light hadrons [25] (see also [26]). It takes into account both 
radiative and collisional energy loss, and fluctuations of the fast parton path lengths in the 
QGP. The calculations of radiative and collisional energy loss are performed with running 
coupling. In our recent note [27] we analyzed within this approach the flavor dependence 
of Raa at the LHC energy. In this paper we concentrate on suppression of single electrons 
and analyze both the RHIC and LHC data. Besides the nuclear modiflcation factor Raa 
we present results for the azimuthal asymmetry V2. 

2 Main features of the model 

In this section we present the basic features of our approach. We refer the interested 
reader to Refs. [25, 26] for more details. 

We deflne the nuclear modification factor Raa for a given impact parameter b as 

„ „ , dN(A + A^ h + X)/dpTdy 

RAA{b,PT,y)- 



TAA{h)da{N + N ^h + X)/dpTdy ' ' 

where pr is the particle transverse momentum, y is rapidity (we consider the central 



region around y = 0), T4a(&) = j dpTA{p)TA{p — b), Ta is the nucleus profile function. 
The numerator in (1) is the differential yield of the process A + A ^ h + X (for clarity 
we omit the argument h) given by 

dN{A + A-^h + X) r dam{N + N^h + X) 

— = / dpTA{p)TA{p - b) — , (2) 

where dam/dpTdy is the medium-modified hard cross section. Similarly to the ordinary 
pQCD formula we write it as 

dam{N + N^h + X) f^dz^^, ^.da{N + N^i + X) 






dpxdy ^ Ja z^ ' dp^rpdy 

where p^ = Pt/z is the transverse momentum of the initial hard parton, 
da{N + N ^ i + X)/dp\^dy is the ordinary hard cross section, D^,- is the medium- 
modified fragmentation function (FF) for transition of a parton i into the observed particle 
h. For the initial virtuality Q we use the parton momentum p^. As in [25], the hard cross 
sections on the right-hand side of (3) were calculated using the LO pQCD formula with 
the CTEQ6 [28] parton distribution functions. The higher order effects were simulated 
taking for the virtuality scale in ag the value cQ with c = 0.265 as in the PYTHIA event 
generator [29]. This prescription gives a fairly good description of the p^-dependence 
of the spectra in pp-collisions ^. We account for the nuclear modification of the parton 
densities with the EKS98 correction [31] (which gives a small deviation of Raa from unity 
even without parton energy loss). 

The formation length arguments allow, in first approximation, to neglect the overlap 
between the DGLAP and the induced stages of the parton showering [25]. Then, assuming 
that formation of the final particle h occurs outside the medium, symbolically the medium- 
modified FF can be written as 

Z^^TaIQ) ^ A/,(Qo) ® Df/, ® D,/,iQ) , (4) 

where ® denotes ^-convolution, D^/i is the ordinary DGLAP FF for i ^ k parton tran- 
sition, DJ]/. is the FF for j ^ k parton transition in the QGP due to induced gluon 
emission, and -D/i/j describes fragmentation of the parton j into the detected particle h 
outside of the QGP. 

We computed the DGLAP FFs using the PYTHIA event generator [29]. For the 
stage outside the QGP for light partons we use for the Dh/j{Qo) the KKP [32] FFs with 
Qo = 2 GeV. We treat the formation of single electrons from heavy quarks as the two-step 
fragmentations c ^ D ^ e and b ^ B ^ e. For the c ^ D and b ^ B transitions we 
use the Peterson FF with parameters ec = 0.06 and Cf, = 0.006. The ^-distribution for 
the M — )■ e transitions (for M = B/D) that we need may be expressed via the electron 
momentum spectrum dB /dp in the heavy meson rest frame as 

p. _ P f^ , 2 cosh(0-g) dB 
A Jo p^ cosh dp 



^For the heavy quark cross sections our LO formulas give the px-dependences (and the cjb ratio) that 
agree weU with the more sophisticated FONLL calculations [30]. However, the normalization of the cross 
sections are smaller by a factor ^1/2. But for Raa it is not important. 



where 

P = J {q^ + ml) cos\? {(p — 6) — ml^ 9 = arcsinh(P/M), = a,TCsmh.{zP/Jq'^ + ml.), P 
is the heavy meson momentum, and M is its mass. We evaluated the B/D — )■ e FFs 
using the CLEO data [33, 34] on the electron spectra in the i?/D-meson decays. We did 
not include the -B — )■ -D — )■ e process, which gives a negligible contribution [30]. 

To calculate the FFs DJ/^ in the induced stage we use the one gluon spectrum dP/dx 
computed in the LCPI formalism [7] with the help of the method suggested in [35]. 
The formulas for calculation of dP/dx are recorded in the Appendix for the reader's 
convenience. The effect of multiple gluon emission was accounted for using Landau's 
method as in [36] (see [25] for details). As in [25], for the quasiparticle masses of light 
quarks and gluon we take m,g = 300 and m,g = 400 MeV supported by the analysis of the 
lattice data [23] . The results are practically insensitive to the light quark mass. For heavy 
quarks we take m,c = 1.2 GeV and rrif, = 4.75 GeV. We use the Debye mass obtained in the 
lattice calculations [37] giving fj,D/T slowly decreasing with T {hd/T ~ 3 at T ~ 1.5Tc, 
Hd/T ~ 2.4 at T ~ 4Tc). We used running Ug frozen at some value q;{^ at low momenta 
(the technical details for incorporating the running as can be found in [35]). For gluon 
emission in vacuum a reasonable choice is af*" ~ 0.7 [38, 39]. However, in the QGP the 
thermal effects can suppress the q;{^, and we regard it as a free parameter which should 
be fixed by the data. If our model is valid, the af'' for light hadrons and single electrons 
should be close to each other. 

As in [25], we treat the coUisional mechanism as a perturbation to the radiative one. 
We account for its effect by redefining the initial QGP temperature in calculating the 
radiative medium-modified FFs according to the condition 

^EradiT^) = AEradiTo) + AE^^l (Tq) , (6) 

where AErad/coi is the radiative/collisional energy loss, Tq is the real initial temperature 
of the QGP, and Tq is the renormalized temperature. We carry out this temperature 
renormalization for each parton trajectory in the QGP (separately for quarks and gluons). 
For the collisional energy loss we use the Bjorken method [11] with an accurate treatment 
of kinematics of the 2 — )■ 2 processes (the details can be found in [22]). For the collisional 
mechanism we use the same parametrization of as{Q) as for the radiative one. Both the 
radiative and collisional contributions in (6) were calculated for maximum energy transfer 
constrained by half of the initial parton energy. 

3 Numerical results and comparison with the data 

We have performed the computations using Bjorken's 1+lD expansion of the QGP [40], 
which gives T^tq = T^t. We take Tq = 0.5 fm. In calculating the medium-modified 
FFs, for simplicity, we neglect variation of the initial temperature Tq with the transverse 
coordinates across the overlapping area of two colliding nuclei. We define this area as 
overlapping of two circles with radius R = Ra + kd, where Ra and d are the parameters 
of the Woods-Saxon nuclear density pAij') = Po/[l + exp((r — _R^)/(i)]. We take k = 1.5, 
which guarantees that the fraction of the lost QGP volume is negligible. The results are 
not very sensitive to variation of k in the physically reasonable range 1 < A; < 2 (except 
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Figure 1: The radiative (solid) and collisional (dashed) energy loss for (a) -u, rf, (b) c, and 
(c) h quarks in expanding plasma of size L = 5 f m for Tq = 300, 400, and 500 MeV for 
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0.5. The order of the curves corresponds to ordering of their Tq. 



for very peripheral AA-collisions, which we do not address in the present paper). To fix 
To (in each centrality bin) we use data on the charged hadron multiplicity pseudorapidity 
density dNch/drj from RHIC [41] and LHC [42, 43]. For the entropy/multiplicity ratio we 
use dS/dy dNch/dr] ^ 7.67 obtained in [44]. For the chemically equilibrated ideal QGP 
(we take Nf = 2.5) it gives Tq ~ 320 MeV for central Au+Au collisions at y/s = 200 
GeV, and Tq ~ 420 MeV for central Pb+Pb collisions at ^/s = 2.76 TeV. For each hard 
parton we calculate accurately the path length in the QGP, L, according to the geometry 
of the AA-collision. Since the QGP should cool quickly at times about 1 -r- 2 units of the 
nucleus radius due to transverse expansion [40], we impose the condition L < Lmax- We 
performed the computations for L^ax = 8. We checked that the bigger value Lmax = 10 
fm gives almost the same. 

In order to illustrate the relative contribution of the collisional mechanism to parton 
energy loss in our model, in Fig. 1 we present AErad and AEcoi for light and heavy quarks 
at To =300, 400 and 500 MeV obtained for L = 5 fm, which is a typical parton path 
length for central AA-collisions at RHIC and LHC. One can see that for light quarks and 
c-quark the contribution of the collisional mechanism is relatively small at p > 5 GeV. But 
for 6-quark the collisional mechanism becomes clearly important at p < 10 GeV. In this 
region our prescription (6) does not apply. From the variation of AE for the radiative and 
collisional mechanisms with the initial plasma temperature in Fig. 1 one can understand 
the magnitude of the ratio Tq/Tq at Tq ~ 300 H- 400 MeV relevant to RHIC and LHC. For 
light quarks and c quark at p ~ 10 ^ 50 GeV (Tq/Tq)^ ~ 1.3 -^ 1.5, and for b quark at 
j9 ~ 10 H- 20 GeV this ratio is somewhat larger ~ 1.5 -^ 1.8. 

Fig. 2 shows comparison of our predictions for Raa for 0-5% centrality bin for (a) 
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Figure 2: (a) Raa for 7r° for 0-5% central Au+Au collisions at ^i = 200 GeV from our 
calculations compared to data from PHENIX [45]. (b,c) Raa for charged hadrons for 0-5% 
central Pb+Pb collisions at ^/s = 2.76 GeV from our calculations compared to data from 
(b) ALICE [46] and (c) CMS [47]. Systematic experimental errors are shown as shaded 
areas. The curves show our calculations for radiative and collisional energy loss (solid), 
and for purely radiative energy loss (dashed) for a{^ = 0.4 (upper curves) and 0.5 (lower 
curves) . 



TT^-meson in Au+Au collisions at ^/s = 200 GeV to PHENIX data [45], and for (b,c) 
charged hadrons in Pb+Pb collisions at y/s = 2.76 TeV to (b) ALICE [46] and (c) CMS 
[47] data. Note that in our calculations based on (4) we ignore possible anomalous baryon 
contribution [48] to the yield of charged particles. Theoretically it is expected to be small 
at LHC [48] . In our model with the KKP FFs [32] for the hadronization outside the QGP, 
Raa for charged hadrons turns out to be very close to that for pions. Experimentally the 
preliminary data from ALICE [49] on Raa for neutral pions also corroborate this. We 
present our results for a^^ = 0.4 (upper curves) and 0.5 (lower curves). To illustrate 
the effect of collisional energy loss we show the total Raa with radiative and collisional 
energy loss (solid) and for purely radiative energy loss (dashed). One can see that the 
effect of the collisional mechanism is relatively small (especially for LHC). We present 
the results for pt >, 5 GeV since for smaller momenta our calculations of the induced 
gluon emission (based on the relativistic approximation) are hardly robust. Fig. 2 shows 
that for light hadrons the window a^"^ ~ 0.4 -^ 0.5 leads to a reasonable magnitude of 
Raa- However, the agreement in the shape of the p^-dependence of Raa is evidently not 
perfect. But this discrepancy does not seem to be very dramatic since the theoretical 
uncertainties of the approximations involved may be significant. One of the most serious 
sources of the theoretical errors, that can be important for the p^-dependence of Raa, is 
the Landau approximation for multiple gluon emission [36]. As can be seen from Fig. 2a 
for RHIC the agreement of the theoretical Raa (radiative plus collisional energy loss) 
with the data is better for a^^ = 0.5. Figs. 2b, c show that for LHC the value a{^ = 0.4 
seems to be preferred by the data (if one considers the complete pt range). Thus, the 
values ai^ = 0.5 and 0.4 seem to be reasonable benchmarks for calculations of the nuclear 
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Figure 3: The electron Raa in Au+Au collisions at y/s = 200 GeV for (a) 0-5%, (b) 10- 
40%, (c) 0-10%, (d) 20-40% centrality classes. The curves show calculations for radiative 
and collisional energy loss (solid), and for purely radiative energy loss (dashed) including 
charm and bottom contributions for a{^ = 0.4 (upper curves) and 0.5 (lower curves). 
Data points are from STAR [2] and PHENIX [3]. Systematic errors are shown as shaded 
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modification factor of single electrons at RHIC and LHC energies. The tendency of the 
decrease of a{^ from RHIC to LHC, first observed in [26], is natural, since the thermal 
reduction of ag should be stronger at the LHC energies. In fact, the variation of a{^ may 
be stronger if one takes Tq^^'" > Tq^'", which seems to be quite reasonable since from 
the dimension arguments one can expect tq oc 1/Tq. However, our purpose is to study 
the variation of nuclear suppression from light flavors to heavy ones probed via single 
electrons, and for this reason it is sufficient to have just a(^ fixed at each energy from the 
light hadron data. 

In Fig. 3 we compare results of our model with STAR [2] and PHENIX [3] data on 
the electron Raa- Comparison to the data from ALICE [4] is shown in Fig. 4a. In 
Fig. 3, 4a we show the total (charm plus bottom) Raa with (solid) and without (dashed) 
collisional energy loss. From Figs. 3, 4a one sees that our pQCD model for the same 
window of a{^ as for light hadrons leads to quite satisfactory agreement with data on the 
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Figure 4: (a) The electron Raa for 0-10% central Pb+Pb collisions at y^ = 2.76 TeV 
for «{'' = 0.4 (upper curves) and 0.5 (lower curves) including charm and bottom. The 
experiment points are the preliminary ALICE data [4]. Systematic errors are shown as 
shaded areas, (b) The electron Raa for charm (thick) and bottom (thin) contributions at 
0.4. In (a,b) the solid curves show results for radiative plus collisional energy loss, 
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and the dashed curves for purely radiative energy loss. 



electron Raa- Similarly to Raa for light hadrons the electron data support af*" ~ 0.5 
for RHIC, and a^'' k. 0.4 for LHC. Thus, the simultaneous description of the nuclear 
suppression of light hadrons and single electrons in the pQCD picture seems quite possible. 
Of course, this should be taken with a caution since the overlapping of the pt region, where 
our approximations make sense, with that studied experimentally is still rather narrow 
(especially for RHIC), and namely in this region the experimental errors are very large. 

To illustrate the effect of collisional energy loss for c — )■ e and 6 — )■ e processes in Fig. 4b 
we show the curves for Raa separately for charm (thick) and bottom (thin) with (solid) 
and without (dashed) collisional energy loss for a{^ = 0.4. One sees that the collisional 
mechanism is more important for bottom. The effect becomes especially significant at low 
Pt- At Pt ^ 5 + 6 GeV our treatment of the collisional mechanism as a perturbation to the 
radiative one, with the help of (6), loses accuracy. Evidently, in this regime the radiative 
and collisional mechanisms must be treated on an equal footing. However, a solution of 
this challenging problem is still lacking. For charm the situation is better, since across 
the whole energy range, where the relativistic approximation makes sense, the collisional 
energy loss remains relatively small. Note that, as one can see from Fig. 4b, the difference 
between the charm and bottom suppression factors at p^ ^ 6 GeV is relatively small. For 
this reason the total (charm plus bottom) Raa is quite stable against variation of the h/c 
ratio which is not very robust. Even at pj- ~ 6 + 10 GeV variation in h/c ratio of ±30% 
changes the total Raa by ±(2 + 3)%. 

It is interesting to compare with the data the theoretical predictions for the electron 
azimuthal anisotropy V2-, which is sensitive to the L-dependence of the heavy quark energy 
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Figure 5: V2 for single electrons for 20-40% centrality class in (a) Au+Au collisions at 
y/s = 200 GeV, and in (b) Pb+Pb collisions at y/s = 2.76 TeV. The theoretical curves 
are for a{^ = 0.4 (lower curves), and 0.5 (upper curves), and include both the charm 
and bottom contributions. Data points are from (a) PHENIX [3], and (b) ALICE [4]. 
Systematic errors are shown as shaded areas. 



loss. In Fig. 5 we compare our calculations for the 20-40% centrality bin to V2 from (a) 
PHENIX [3] and (b) ALICE [4]. One sees that the agreement with the ALICE data is 
fairly good. However, the experimental errors are very large and the px range is too 
limited to make a definitive conclusion on the preferred value of «{''. For the PHENIX 
data [3] Pt ^ 4 GeV. For such low pt our calculations are not robust. Nevertheless, our 
V2 for a{.^ = 0.5 (favored by data on R^a for pions) at pt = 6 GeV matches reasonably 
well the experimental f 2 at p^ ~ 4 GeV. 

Note that our predictions were obtained with the radiative energy loss for the QGP 
modelled by a system of the static Debye screened color centers [5] . The generalization 
to the dynamical QGP described within the HTL scheme is trivial (see [50] for details). 
It is reduced to replacement of the potential (A4) (let us call it Vgtat) by a dynamical 
potential Vdyn, which can be expressed through the gluon polarization operator. In the 
HTL scheme an elegant formula derived in [51] allows to write Vdyn similarly to the static 
case just replacing the factor l/(g^ + fioY by l/[g^(g^ + fin)] in the formula for the 
dipole cross section (A6), and increasing a little the overall normalization by a factor 
^f^(l + Nf/Q)/{1 + Nf/A) ^ 1.19 (for Nf = 2.5). This modification leads to unlimited 
growth of Vdyn at large p (due to zero magnetic mass in the HTL approximation), while 
Vstat flattens at p > l/mD- It has been recently claimed [18] that the dynamical effects 
enhance the heavy flavor suppression, and are important for description of the single 
electron suppression. Of course, modification of the potential in the above manner for the 
dynamical QGP should enhance the heavy quark suppression. However, this enhancement 
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itself is uninteresting in the context of the heavy-to-light ratio R 

dynamical potential should enhance the light flavor suppression as well. One can expect 
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that the dynamical effects affect the induced gluon emission for hght and heavy flavors 
similarly. Indeed, the typical p, dominating the induced gluon spectrum, are the same for 
light and heavy flavors^, it is the region p < y/Su^/rrig [7] (here Slpm is the suppresion 
factor due to the Landau-Pomeranchuk-Migdal effect, which is typically not very small 
(~ 0.2 -^ 0.5) for RHIC and LHC). In this p-region the shapes of Vdyn{p) and Vstat{p) are 
qualitatively very similar. For this reason modification of the nuclear suppression due 
to the dynamical effects should be similar for light and heavy flavors. And, once the 
coupling is fixed by the data on Raa for light hadrons, the dynamical formulation should 
give Raa for heavy flavors (and single electrons) close to that for the static model. We 
have checked this performing numerical calculations (for a fixed coupling constant as in 
the HTL scheme), and have found a negligible modification of the heavy-to-light ratio for 
the dynamical model. Note that the use of the Vdynip) hardly makes the calculations more 
robust. Indeed, the difference between Vdyn{p) and Vgtatip) at p > l/rriD is related mostly 
to the zero magnetic mass in the HTL scheme. But the lattice calculations [53] show 
that in reality the magnetic mass may be of the order of the electric Debye mass. On the 
other hand, in the region p ^ l/m£), where the HTL approximation is not supposed to 
be valid, it gives incorrect normalization of the potential (contrary to the static model, 
which gives the correct result at p — )■ 0). 

4 Summary 

In this paper we have examined whether it is possible in the pQCD picture of the parton 
energy loss to describe simultaneously experimental data from RHIC and LHC on nuclear 
suppression of light hadrons and single electrons from the heavy meson decays. We 
have performed calculations taking into account radiative and collisional energy loss, and 
fluctuations of the fast parton path lengths in the QGP. The calculations of radiative 
energy loss have been performed within the LCPI approach [7]. Both radiative and 
collisional energy loss were computed with running as frozen at low momenta at some 
value a{^. 

We have found that once the value of a{^ is fixed from the data on the nuclear mod- 
ification factor for light hadrons it gives a satisfactory agreement with the data on the 
electron Raa as well. For the electron azimuthal anisotropy V2 the agreement is also 
within the experimental errors. Our calculations show that the effect of collisional energy 
loss is realtively small, and cannot be crucial for the flavor dependence of the nuclear 
suppression factor for pt > 10 GeV. The collisional mechanism becomes very important 
only for the bottom contribution to the electron spectrum at momenta < 6 -i- 8 GeV. 

Our results, together with fairly good agreement with the ALICE LHC data on Raa 
of D-mesons [54, 55] obtained in our recent analysis [27], give support for the pQCD 
picture of parton energy loss both for light and heavy flavors. However, the available 
data on the nuclear suppression of single electrons, especially from RHIC, analyzed in the 
present work, are restricted to rather low px, where the conditions of applicability of our 



^For heavy quarks the typical p becomes smaller than for light partons only in the tail region x > 
mg/mQ, where quark mass suppresses strongly gluon emission. 
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pQCD model may be not good enough. For more conclusive test of the model it is highly 
desirable to have data on the electron Raa that extend to larger values oipT- 
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Appendix: Formulas for one gluon x-spectrum 

We use the representation of the one gluon emission ^-distribution obtained in [35] which 
is convenient for numerical calculations. For q ^ gq process it reads 



dP 
dx 



dzn{z) 



daf/}{x,z) 



dx 



(Al) 



where n{z) is the medium number density, da^/i/dx is an effective Bethe-Heitler cross 
section accounting for both the LPM and finite-size effects. The da^f^Jdx reads 



daf/}{x,z) 



dx 



^'n™/.WQ^(0)^^^«-^)' 



ttM 







dp\^ 



(A2) 



p=0 



Here PHx) = Cf[1 + (1 — x)'^]/x is the usual splitting function for q ^ gq process, 
M = Ex{l — x) is the reduced "Schrodinger mass", Q^(0 = ciM/^ with a ^ 1.85 [22], 
F is the solution to the radial Schrodinger equation for the azimuthal quantum number 
m = 1 






1 



d 



2M \dp 



. , Am^ - 1 1 



F{C,p) (A3) 



with the boundary condition F(^ = 0,p) = y/pas{p,x, z)eKi{ep) (Ki is the Bessel func- 
tion), Lf = 2M/e^ with e^ = m'^x'^+myi—xy, a^lp, x, z) is the cross section of interaction 
of the qqg system with a medium constituent located at z. The potential v in (A3) reads 

.n{z)a3{p,x,z) 



v{p,x,z) 



(A4) 



The (T3 is given by [52] 



9 1 

a^ip, X, Z) = -[cTgqip, Z) + Orgg-((l - X)p, Z)\ - -aqq{xp, z) 



(A5) 



where 



(Tqq{p,z) = CtCf / dqa^{q 



2/ 2A[l-exp(^qP)] 



(A6) 



is the local dipole cross section for the color singlet qq pair {Cf,t are the color Casimir 
for the quark and thermal parton (quark or gluon), po is the local Debye mass). 



11 



For g ^)- gg one should replace the splitting function and rriq by rrig in e^. The o"3 in 
this case reads 

9 
as{p, X, z) = -[(Tggip, z) + (yqq{{\ - x)p, z) + agg-(a;p, z)\ . (AT) 
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